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Abstract: 

Coupled  surface  plasmon  modes  of  graphene  on  a  plasma  layer  are  studied.  Three 
coupled  plasmon  modes  are  found,  and  analytical  equations  for  their  dispersion 
characteristics  are  given  in  the  long-wavelength  limit.  Symmetric  and  antisymmetric 
modes  are  identified;  between  two  strongly  coupled  modes,  the  antisymmetric  mode  is 
found  to  be  more  energetic  than  the  symmetric  mode.  The  great  tunability  of  plasmon 
characteristics  is  shown  by  changing  the  surrounding  dielectric  constant,  the  coupling 
distance,  and  the  plasma  frequency  of  the  substrate.  (Appl.  Phys.  Lett.  103,  201104 
(2013)) 

Based  on  the  understanding  of  coupled  surface  plasmon  modes  of  graphene,  a 
graphene-based  terahertz  plasmonic  waveguide  is  proposed.  The  proposed  structure 
benefits  from  the  enhanced  confinement  and  increased  attenuation  length  of  graphene 
surface  plasmon  by  placing  the  graphene  sheet  in  proximity  of  metal  layers.  (Appl.  Phys. 
Lett.  105,  011604  (2014)).  A  graphene  plasmon-based  THz  device  is  also  fabricated  and 
measured.  The  experimental  results  show  a  consistent  trend  between  the  plasmon 
resonance  frequency  and  the  grating  period,  as  predicted  by  our  theoretical  analysis. 
(Submitted  to  Sci.  Rep.  (2016)) 


I  ntroduction: 

Graphene  is  a  one-atom-thick  2D  system  that  has  a  unique  hexagonal  crystal  structure 
of  two  carbon  atoms  per  unit  cell.  Unlike  any  other  2D  semiconductor  material  known 
today,  intrinsic  graphene  has  a  zero  bandgap  with  its  charged  carriers  behaving  like  Dirac 
fermions  with  a  zero  mass,  resulting  in  many  extraordinary  properties  that  are  very 
different  from  other  materials.  Due  to  the  lack  of  the  understanding  of  the  graphene 
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plasmons  and  their  interaction  with  the  immediate  surroundings,  we  investigated  the 
coupled  surface  plasmon  modes  of  graphene  and  graphene-based  plasmonic  devices. 

Recently,  we  derived  analytical  models  using  Maxwell  equations  for  coupled  plasmon 
modes  of  structures  that  consist  of  a  sheet  of  graphene  and  a  plasma  substrate  such  as  a 
metal,  a  doped  semiconductor,  or  another  graphene  layer.  Using  these  models,  various 
plasmonic  devices  with  desirable  characteristics  can  be  designed,  such  as  graphene-based 
waveguides  and  modulators.  Examples  of  graphene-metamaterial  hybrid  system  is  shown 
in  Fig.  1,  where  e  ,  e.m,  e,,  e7,  and  e7  are  the  permittivities  of  the  dielectric  gap  filling, 

the  medium  that  is  used  to  construct  the  metamaterial,  the  medium  below  the 
metamaterial,  the  medium  below  graphene,  and  the  medium  above  graphene,  respectively; 
dx  is  the  thickness  of  the  metamaterial;  d?  is  the  distance  between  the  graphene  sheet 

and  the  metamaterial;  other  structural  parameters  charactering  the  miniature  patterns  of 
the  metamaterial  are  also  specified,  such  as  the  width  w  and  the  height  h  of  the  engraved 
grooves.  Depending  on  the  design  of  the  device,  different  types  of  the  surface  excitations 
can  be  induced.  These  excitations  in  turn  give  the  desired  functionality  of  the  devices. 


Fig.  1.  Various  excitations  of  graphene-based  THz  devices. 


Experiment,  Results  and  Discussion: 

(1)  Coupled  surface  plasmon  modes  of  graphene  in  close  proximity  to  a  plasma  layer 
(Appl.  Phys.  Lett,  103,  201104  (2013)) 

It  has  been  shown  that  it  is  possible  to  fundamentally  change  the  plasmon  dispersion 
curve  if  substrate-induced  effects  are  considered  or  if  graphene  is  in  close  proximity  with 
another  plasma  material,  such  as  another  layer  of  graphene,  a  metal  layer,  or  a  metal 
particle.  The  characteristic  plasmon  dispersion  for  the  coupled  plasmon  modes  of 
graphene  that  is  30  nm  away  from  a  doped  semiconductor  thin  film  are  shown  in  Fig.  2. 
Three  coupled  plasmon  modes  &>spl ,  cosp2  and  ojg  are  found  with  to  >  >  rog  (thick 

solid  curves).  The  decoupled  modes  are  also  plotted  as  dashed  curves  for  dl=d2=5  pm 
between  different  layers  of  the  structure;  these  modes  are  decoupled  as  long  as  the 
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plasmonic  wavelength  ql  d...  Two  dashed  lines  are  the  surface  plasmon  modes 

associated  with  the  top  and  bottom  surfaces  of  the  semiconductor  substrate.  Another 
rising  dashed  curve  is  the  graphene  plasmon  mode.  When  q\  surface  plasmon 

fields  on  the  graphene  sheet  and  on  the  semiconductor  surfaces  can  hardly  “see”  each 
other;  in  this  situation,  the  plasmon  modes  are  weakly  coupled,  following  the  dashed 
curves  in  the  high-g  region.  By  contrast,  the  plasmon  modes  are  strongly  coupled  in  the 
limit  q  — >  0  when  retardation  is  neglected,  resulting  in  the  solid  dispersion  curves  that 
significantly  deviate  from  the  dashed  curves.  This  situation  is  similar  to  the  case  of  planar 
metallic  waveguides  where  the  surface  plasmon  mode  on  each  surface  are  strongly 
coupled  when  q  — >  0. 


Fig.  2.  (a)  Coupled  plasmon  modes  (<»  ,<»  and  solid  curves)  of  graphene  30  nm 
away  from  a  semiconductor  thin  film  of  50  nm  thickness,  (b)  Coupled  plasmon  modes 
( a>  and  solid  curves)  of  graphene  30  nm  away  from  a  semi-infinite  semiconductor 

substrate  ( d^  — >  qo).  Decoupled  modes,  calculated  with  d2=  5 pm ,  are  also  plotted  as 
dashed  curves. 

(2)  Enhanced  graphene  plasmon  waveguiding  in  a  layered  graphene-metal  structure 
(Appl.  Phys.  Lett,  105,  011604  (2014)) 

In  our  previous  work,  we  have  shown  that  the  surface  plasmon  dispersion  of  graphene 
is  strongly  modified  by  the  presence  of  a  plasma  material  such  as  a  metal.  In  this  study, 
we  further  show  that  within  a  certain  spectral  region  the  losses  can  be  reduced  when 
graphene  is  placed  near  a  metal  surface.  The  dispersion  characteristics  of  surface 
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plasmons  for  the  structure  shown  in  Fig.  3(a)  are  plotted  in  Fig.  3(c).  The  curves  are 
calculated  using  different  values  of  d2=d3  =  d  marked  next  to  each  curve.  As  can  be 
seen,  the  approximate  analytical  results  (dashed  curves)  agree  very  well  with  the  exact 
numerical  results  in  the  low-frequency  region  as  long  as  the  condition  qxd  1  is 
satisfied.  In  this  region,  the  plasmon  dispersion  is  a  strong  function  of  d;  at  a  fix 
frequency,  ql  increases  with  decreasing  d,  signifying  that  the  confinement  of  the 

plasmonic  field  can  be  consistently  enhanced  with  decreasing  graphene-metal  distance. 
However,  in  the  high-frequency  region,  all  numerical  curves  merge  into  the  blue  curve 
because  of  the  strong  confinement  of  the  plasmonic  field  on  the  graphene  sheet:  the 
presence  of  metal  slabs  becomes  insignificant  as  the  plasmonic  wavelength  becomes 
much  smaller  than  the  physical  distances  between  the  graphene  sheet  and  the  metal 
surfaces. 

The  attenuation  length  q~ 1  is  plotted  in  Fig.  3(d)  with  the  graphene-metal  distance  d 

marked  next  to  each  curve.  As  can  be  seen,  the  approximate  analytical  results  (dashed 
curves)  agree  very  well  with  the  exact  numerical  results  in  the  low-frequency  region.  In 
the  intermediate-frequency  region  where  q~x  is  comparable  to  d,  the  black  curves  are 

below  the  blue  curve,  signifying  reduced  loss  in  the  presence  of  the  metal  slabs. 
Therefore,  as  seen  in  Fig.  3(c)  and  Fig.  3(d),  we  find  that  an  enhanced  confinement  and 
an  increased  propagation  distance  can  be  simultaneously  achieved  by  the  presence  of  the 
metal  slabs  within  a  certain  spectral  region. 


(a)  (b) 
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Fig.  3.  (a)  Planar  waveguide,  and  (b)  nonplanar  rectangular  graphene  waveguide,  (c) 
Dispersion  characteristics  and  (d)  the  attenuation  length  of  surface  plasmons  as  a  function 
of  surface  plasmon  energy  for  the  structure  shown  in  (a)  with  d  =d  =  d. 
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(3)  Family  of  graphene-assisted  resonant  surface  optical  excitations  for  terahertz 
devices  (submitted  to  Sci.  Rep,  (2016)) 


Frequency  (THz)  f  Frequency  (THz) 


Fig.  4.  (a)  Images  of  scanning  electron  microscope  for  ID  and  2D  gratings,  (b)  Ratio  of 
the  real  part  of  the  transmitted  field  through  gratings  to  that  of  the  free  space  without  the 
grating,  (c)  and  (d)  are  the  simulated  and  experimental  data  of  the  transmitted  field 
through  graphene-grating  structure,  respectively,  (e)  Transmitted  fields  for  different 
values  of  the  polarization  angle,  (f)  Transmittance  measurement  of  FTIR  for  ID  gratings 
with  (thick  curves)  and  without  (thin  curves)  graphene  on  top. 


The  majority  of  the  proposed  graphene-based  THz  devices  consist  of  a  metamaterial 
that  can  optically  interact  with  graphene.  This  coupled  graphene-metamaterial  system 
gives  rise  to  a  family  of  resonant  modes  such  as  the  surface  plasmon  polariton  (SPP) 
modes  of  graphene.  In  this  study,  graphene  SPP-based  THz  devices  are  fabricated  and 
measured,  as  shown  in  Fig.  4(a).  Terahertz  time-domain  spectroscopy  (THz-TDS)  is  used 
to  measure  the  transmitted  electric  field  through  the  grating.  The  experimental  and 
simulation  results  are  plotted  in  Fig.  4(c)  and  Fig.  4(d),  respectively.  As  can  be  seen,  the 
plasmonic  behavior  is  well  described  by  the  simulation.  By  rotating  the  angle  from  0°  to 
90°  for  the  2D  sample,  as  plotted  in  Fig  .4(e),  a  smooth  transition  is  shown  for  the  curves 
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at  different  angles,  indicating  an  angle-dependent  scattering  rate  and  the  polarization- 
sensitive  nature  of  our  2D  structure.  The  effect  of  the  Drude  conductivity  is  observed 
using  Fourier-transform  infrared  spectroscopy  (FTIR),  shown  as  the  shaded  areas  in  Fig. 
4(f).  Because  the  graphene  SPP  mode  is  supported  by  the  Drude  conductivity,  the 
operational  range  of  the  graphene  SPP-based  THz  device  shown  is  limited  not  only  by 
intrinsic  Landau  damping,  but  also  by  the  frequency-dependent  nature  of  the  Drude 
conductivity  of  graphene. 
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